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Abstract My personal views of the Chiral Magnetic Effect are presented, which 
starts with a story about how we came up with the electric-current formula and 
continues to unsettled subtleties in the formula. There are desirable features in the 
^ C| formula of the Chiral Magnetic Effect but some considerations would lead us to 

Q_[ even more questions than elucidations. The interpretation of the produced current 

is indeed very non-trivial and it involves a lot of confusions that have not been 
resolved. 

1 Introduction - Discovery of the Chiral Magnetic Effect 

The Chiral Magnetic Effect (CME) is concisely summarized in the following handy 
formula; 

which represents an electric current associated with the non-zero chirality and the 
external magnetic field B. Here stands for the number of colors in quantum chro- 
modynamics (QCD) and qf represents the electric charge carried by the quark flavor 
/ where / runs over up, down, strange, etc. Equation (Q~|) looks simple, but the phys- 
ical meaning of this CME current is far from simple. Let me begin with telling some 
historical remarks on the discovery of the CME-current formula, hoping that it may 
be instructive and even inspiring to some readers. 

When we, Harmen Warringa, Dima Kharzeev, and I, started working on the com- 
putation of j, we had no a priori idea about the final answer, hence we did not really 
expect that the final result should be such elegant and beautiful. For several years 
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Harmen and Dima had been working on the implication of axion physics in the con- 
text of the relativistic heavy-ion collision [I will come back to the relevance of the 
CME to axion physics later] At that time, around the year of 2007, I was think- 
ing of a different (but related) physics problem, i.e. color-superconducting states 
in a strong B inspired by a pioneering work |8|. Harmen and I just chatted in the 
corridor of the RIKEN BNL Research Center (RBRC) about B-effects on color su- 
perconductivity, which was soon promoted to intriguing discussions, and a fruitful 
collaboration after all. Nearly simultaneously with the successful completion of our 
project on color superconductivity in B (16 1 (see also Ref. [34] for an accidental 
coincidence of the research interest with our Ref. ifTBI ). a monumental paper by 
Harmen, Dima, and Larry McLerran appeared ll28l . While we were finalizing the 
color-superconductivity paper (or struggling with referees, probably), Harmen ex- 
citedly explained the idea of the Chiral Magnetic Effect to me. Also, I clearly re- 
member that Larry came over mischievously (as always) to ask about the strength of 
my B in the neutron-star environment (eB ~ 10 15 gauss at most on the magnetar sur- 
face). As compared to their B produced in the relativistic heavy-ion collision where 
eB ~ 10 20 gauss could be reached, mine was only negligible... Indeed, historically 
speaking, the recognition of such B as strong as the QCD energy scale Aqcd m real- 
istic circumstances was an important turning point to get the B-physics research into 
gear. In other words, physics researches at eB ~ Aq CD have come to make pragmatic 
sense rather than purely academic one since this turning point in 2007. There was 
really a tremendous change in the attitude of researchers. 

One year later, Harmen invited me to his continued project with Dima on the 
Chiral Magnetic Effect. In their first paper the formula was given in a different style 
from what is known today, namely, it was not the current but the charge separation 
Q expressed as [28| 

G = 2e„£|«/|7(2|?/*l). (2) 

/ 

Here Q w is the topological charge (i.e. counter part of ^5 in Eq. (fl~|l) and y{x) 
is a function dependent on the microscopic dynamics of quark matter. Accord- 
ing to the analysis in Ref. 1281 one can approximate y(x) by a simple function; 
y(x < 1 ) = x, y(x > 1 ) = 1 ■ This means that, if the magnetic flux per unit topologi- 
cal domain, <t>, is large enough, Q sa 2Q w Y,f \lf\ - This result is naturally understood 
from the index theorem, i.e. 2Q W = N5 = Nl—Nr. Under such strong B, all the spin 
directions should completely align in parallel with B, and thus the momentum di- 
rections are uniquely determined in accord with the chirality. All produced chirality 
should contribute to the charge separation, leading to Q ~ NsY,/ \<lf\ that is nothing 
but 2(2vf L/ \<lf |- In the weak field case, on the other hand, Q ks 44><2 V vL/9/ was the 
theoretical estimate. 

Equation (f2| is as a meaningful formula as Eq. (HJ, but the determination of y(x) 
requires some assumptions. Besides, since the formula involves Q w , it is unavoid- 
able to think of topologically non-trivial gauge configurations. As a matter of fact, 
Harmen and I once tried to compute Q concretely on top of the real-time topologi- 
cal configuration, namely, the Liischer-Schechter classical solution 13111401 . which 
turned out to be too complicated to be of any practical use. Then, Harmen hit on 
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a brilliant idea to deal with Q w , or strictly speaking, an idea to skirt around Q w . 
[He invented another nice trick later to treat Q w more directly. I will come to this 
point later.] The crucial point is the following; it is not the topological charge Q w 
but the chirality N5 that causes the charge separation. It is tough to think of Q w , then 
what about starting with N5 not caring too much about its microscopic origin? If 
one wants to fix a value of some number, one should introduce a chemical potential 
conjugate to the number. In this case of N5, the necessary ingredient is the chiral 
chemical potential /I5 that couples the chiral-charge operator i//y°y 5 i//. In my opin- 
ion the introduction of /I5 was a simple and great step to make the CME transparent 
to everybody. In this way the CME has eventually gotten equipped with enough 
simplicity and clarity. 

The remaining task was to answer the following question; what is j in a system 
with both /I5 and B? Harmen and I were first going to calculate the expectation 
value of the current operator xfry^W directly (see the derivation A in Ref. ITTlD . To 
this end we had to solve the Dirac equation in the presence of fis and B to construct 
the propagator. Now I am very familiar with the way how to do this explicitly, but 
when we started working on this project, we had not had enough expertise yet, apart 
from some straightforward calculations in color superconductivity. Some years later 
Harmen, Dima, and I wrote a paper in which we reported the diagrammatic method 
to derive Eq. ((TJ (see Appendix A in Ref. [ 12 1). Let me briefly explain this derivation 
here; the electric current in the z-axis direction is written in terms of the propagator 
as 



/ 



xtr 



fPn M 7T^ Pn M 



(3) 



where the /^-integration is either at T = or the Matsubara sum at T ^ 0. If we 
choose the gauge as Aq = A x =A Z = and A y = Bx, the tilde momentum in the de- 
nominator is p = (po,Q, —sgn(qB)y/2\qB\n,p z ). We do not need the explicit form 
of the Landau wave-functions P„ (x) that take a 4 x 4 matrix structure in Dirac space. 
Because we are interested in j || B here, y z commutes with P„(x) and thus we need 
only P n (x) 2 which equals 1 for n > and (1 + isgn(qfB)Y*Y y )/2 f° r n = 0- Af- 
ter some calculations one can confirm that Eq. (O is reduced to Eq. fl} regardless 
of the temperature T and the flavor-dependent mass Mr. Let me make a comment 
on this rather naive calculation. In most cases the proper-time method is the best 
way to proceed in theoretical calculations [42 22] and the above form of the quark 
propagator is not widely known. For the purpose of calculating a finite quantity like 
the CME current, I would like to stress that the above quark propagator should be 
equally useful. Actually it is almost obvious in Eq. (O that any contributions from 
the Landau non-zero modes are vanishing and the current arises from the Landau 
zero-mode only. 

Coming back to the story of our first attempt to discover j, I remember that Har- 
men and I came to the office and brought different answers every morning and had 
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the hottest discussions all the day. It took us a few days until we eventually con- 
vinced ourselves to arrive at the right answer. Later on, Harmen had great efforts 
to dig out several independent derivations of Eq. ([TJ while preparing for our paper. 
Among various derivations we first found the one based on the thermodynamic po- 
tential (i.e. the derivation C in Ref. [11]). Because this calculation plays some role 
in later discussions on the physical interpretation of the CME current, let us take a 
closer look at the detailed derivation using the thermodynamic potential. 

The most essential ingredient is the quasi-particle energy dispersion relation in 
the presence of B and /I5. For B along the z-axis, one can solve the Dirac equation 
to find the dispersion relation, 



CO 



p,s 



(p z z +2\q f B\n)^ 2 + sgn(p z )sn 5 



-M 



f ' 



(4) 



where s is the spin, qj and Mf are the electric charge and the mass of quark flavor 
/. Once the one-particle energy is given, one can immediately write the thermody- 
namic potential down as 



/ .v=± n=0 



dpz 
, ~7m 



co„ 



r£in(i 

± 



-(£Op,s± 



(5) 



at finite temperature T and quark chemical potential fi. The spin factor, a„. s , is 
defined as a„ )S = 1 (n > 0), 8 S+ (n = 0,qfB > 0), <5 4 _ (n = 0,q/B < 0). This factor 
is necessary to take care of the fact that the Landau zero-mode (n = 0) exists for 
one spin state only. The current j z is obtained by differentiating £2 with respect to 
the vector potential A z . Because the vector potential in the matter sector resides only 
through the covariant derivative, the following replacement is possible inside of the 
/^--integration, 

q— ■ (6) 



dA 7 



The combination of this derivative and the ^--integration ends up with the surface 
terms. It is the characteristic feature of the quantum anomaly that a finite answer 
results from the ultraviolet edges in the momentum integration. That is, the CME 
current reads, 



if\if B \ 
2% 



t±n J —A 2% dp z L ^ 



^I^[(A ±M5 )-(A TM5 )]=Ar c £^. 



/ 



/ 



2k 2 



(7) 



Here, in the second and the third lines, ± appears from the Landau zero-mode al- 
lowed by a„ t s, i.e. ± amounts to sgn(qyfi) which cancels the modulus of \q/B\, 



Views of the Chiral Magnetic Effect 



5 



and the matter part drops off for infinitely large (0 lhS {p z = ±A ). It would be just a 
several-line calculation to make sure that Eq. (01 is equivalent with Eq. ^ and they 
are calculations at the one-loop level. It is also a common character of the quantum 
anomaly that the one-loop calculation would often give the full quantum answer. Al- 
though I do not know any explicit check of the higher-order loop effects, the above 
method at the one-loop level is my favorite derivation of Eq. (HJ; all the calculation 
procedures are so elementary and transparent. 



2 Chiral Separation Effect 

Soon later, Harmen and I found that a very similar topological current had been 
discovered in the neutron-star environment, that is, the axial current associated with 
the quark chemical potential fi and the magnetic field B [32], 



This is a chiral dual version of Eq. ((T). Nowadays people call Eq. (0 the Chiral 
Separation Effect (CSE) in contrast to Eq. ([TJ referred to as the Chiral Magnetic 
Effect. When we learned the fact that Eq. ([H} had been known earlier, our excite- 
ment got cooled down a bit. Also, three years later, we came to know that the CME 
formula had been discovered further earlier. Now there is a consensus in the com- 
munity that the CME formula ([T} was first derived by Alex Vilenkin |47l . It was an 
embarrassment for me to have overlooked his work until he brought our attention to 
his old papers. In fact an equivalent of Eq. (HJ has been rediscovered over and over 
again l20l[T9l [T l and I would not be surprised even if Eq. ([TJ is still buried in further 
unknown works. [I am not talking about the recent activities to derive Eq. ([TJ from 
a deeper insight into physics such as Berry's curvature B31 1501 . hydro or kinetic 
approaches [25, 18, 46 24|, and so on, which really deserve more investigations.] 

The derivation of Eq. (O is worth discussing here. The topological effects in 
quantum electrodynamics (QED) from A^ c x Nf quarks add terms in the action as 



associated with an axial rotation by 0(x). In this way we see that the axial current is 
not conserved but anomalous. With the replacement of Aq = ji and e 0,J djA^ = B', 
one can transform this expression using the integration by parts into 




(8) 




(9) 
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from which Eq. ([8]l immediately follows. This derivation also tells us that the B- 
induced current in the right-hand side of Eq. ([8]) is nothing but a part of the Chern- 
Simons current ~ £^ vpa A v d p A c in QED. It should be mentioned that the derivation 
presented above is a little bit cooked up by me for the illustration purpose and one 
should refer to the original paper 1321 for more careful treatments of the surface 
integral. 

Before going on our discussions, let me point out that the above derivation im- 
plicitly assumes massless quarks. If quarks are massive, Eq. (0 should be modified 
with an additional term 2iMf(\jff'y 5 Yf)- This modification would be harmless as 
long as the pseudo-scalar condensate is vanishing, but in principle, Eq. ((HJ could be 
dependent on Mf unlike Eq. ((TJ as argued explicitly in Ref. l32l . In fact it is quite 
subtle whether Eq. (0 is sensitive to Mf or not, and I will address this question in 
an explicit way soon later. 

It would be an interesting question how to derive Eq. dS) microscopically just like 
the ways addressed in the previous section. In fact I have once tried to prove Eq. (H) 
based on the thermodynamic potential by inserting an axial gauge field. There must 
be a way along this line, but I could not solve it (or I would say that I did not have 
enough time to find it out...). Instead, here, let me introduce another derivation based 
on the propagator as in Eq. (|3). 

The axial current is expressed as 



K X g/lg/gl V- 

2n L 



T dpo f dp z 
In J 2n 



dx 



xtr 



ffP n (x)- 



P n {x) 



(11) 



at finite quark chemical potential jx. It is easy to see that any contributions from 
vanish due to the Dirac trace. Only the Landau zero-mode produces a term 
involving y x y y which makes ti('/ ) y c y y y z y s ) = —4/ ^ 0. Then, the above expression 
simplifies as 



^ 2n 



f 



u 2n 



f 



T dpo 
2n 

T dpo 
2n 



2n 



rr {p + H)2-pZ-Mf r 



^Jpo+iif-pl-Mi 



c2 f2n dn 



(12) 



where Z(ju) denotes the partition function at finite density in (l+l)-dimensional 
theory (as a result of the dimensional reduction with the Landau zero-mode), and 
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thus the /i-derivative leads to the quark density n. In the second line we used 
2(P0 + Al)/[(P0 + A0 2 -P?-^] - {d/dn)\n[{p {) + n) 2 -p\-M}]. One might 
have thought that it is a simple exercise to evaluate Z(ji) with the (l+l)-dimensional 
integration. The fact is, however, that the finite-/^ system in (1+1) dimensions is by 
no means simple. 

In Ref. 1321 one can find exactly the same expression as above in a slightly dif- 
ferent calculation and the density is written as (see Eq. (37) in Ref. [32]), 



dp z 
In 



1 



? (w f -ll)/T + j e ((0 f +n)/T + l 



(13) 



with (0/ — J pl+Mj. This result is certainly M^-dependent as suggested in the 
paragraph below Eq. ( fTOb . and this would make a sharp contrast to the CME cur- 
rent ([]])■ 

We know, however, that the density in the (l+l)-dimensional fermionic theory 
arises from the anomaly [41 1 and the density dT3T > is not the right answer. In fact, 
in view of the second line of Eq. (TT2l) . it seems at a glance that the jj. -dependence 
could be absorbed in the /?o-integration, which already gives us an impression that 
something non-natural should be happening. To see this, let us take one-step back 
to the microscopic expression, i.e., the (l+l)-dimensional partition function reads, 

T ^J^tx[f(p + n)-fp z -M f ], (14) 

from which the ji -dependence could be eliminated by the chiral rotation (for the 
zero-mode basis only), 

W^'^Vo, (15) 

leading to (here, we shall show results at T = for simplicity, but nothing is changed 
even at finite T), 

A -v dp, 1 . f A +f* dp, 1 . 
-A+ju 2% 2 J-a-h 2% 2 



with (Of = \jp\ + \Mf\ 2 , where Mf = Mfe 11 "^"^^ is the chirally tilted mass. The 
momentum integration is shifted according to the chiral rotation (TT3T >. The first (sec- 
ond) integral corresponds to the particle (anti-particle, respectively) contribution. 
Thus, one can extract the jj. -dependent piece from the surface terms as follows; 
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Z = 




dx 2 J a 1% 
A* 2 

- — V (/^-independent terms) , 



(17) 



which results in the density n/ = that is independent of Mf 1411 . It is clear 
from the second line of the above calculation that the density originates from the ul- 
traviolet edges, which is the characteristic feature of the anomaly. The full quantum 
answer is then given as 



Equivalently, if one is interested in deriving the same answer from Eq. (fTSl i directly, 
one should split the composite operator as \jf(x)'fy s \lf(x) — > ^{x + e)fy 5 yi{x) and 
insert the infinitesimal gauge connection from x to x + £. Interestingly, contrary to 
Ref. [32], the Chiral Separation Effect (O is presumably insensitive to the quark 
mass just like the Chiral Magnetic Effect (|T). Whether Eq. (JHJ is robust or not re- 
gardless of Mf is an important question particularly in the context of the Chiral 
Magnetic Wave (CMW) [ 30 ] . The anomalous nature of the density (TT~8T > implies 
that the CMW can exist also in the chiral-symmetry breaking phase where quarks 
acquire substantial mass dynamically. 

I would not insist that I could prove the non-renormalization of Eq. <j8j since the 
above is just a one-loop perturbation and non-perturbative interactions may change 
the story; I would like to thank Igor Shovkovy for raising this unanswerable but 
unforgettable question. The interested readers may consult Refs. IF2TI [T5 I for some 
examples of non-non-renormalization. Anyway, I can at least say with confidence 
that, if B is super-strong, the reduction to the (l+l)-dimensional system should be 
strict, and then Eq. ( fT3l must be altered, conceivably as iif — \x.j n l9l. 

Although the interpretation of Eq. <fl~8T > may swerve a bit from our main stream, 
I would emphasize that Eq. < TT~8T > is extremely interesting and it would be defi- 
nitely worth revisiting its profound meaning. Actually, Eq. (Q3) has an impact on 
the structure of the QCD vacuum. Let us consider the hadronic phase with sponta- 
neous breakdown of chiral symmetry. After the rotation 1151 . apart from the anoma- 
lous term ji 2 / (2n), the system is reduced to that at zero density, which means that 
X = (WqWq) should take a finite value. Therefore, in terms of the original fields y/o, 
the chiral condensates form a spiral structure, 



rif = 



Jz y2n 2 



(18) 



(VWo) = Zcos(2jiiz) , (yaffvo) = %sin(2/iz) , 



(19) 



which is called the chiral spiral or the dual chiral-density wave (if y 5 is involved) G] 
[33l . In particular, if the above type of the inhomogeneous ground state is caused by 
B, it is sometimes called the chiral magnetic spiral [5 1. 
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I should emphasize that Eq. ( fT8l ) does not really require the dimensional reduc- 
tion, while the chiral magnetic spiral needs the pseudo (l+l)-dimensional nature 
under sufficiently strong B. This point might be a bit puzzling. As long as jf is con- 
cerned, only the Landau zero-mode remains non-vanishingly for any B, and the mo- 
mentum integration is purely (l+l)-dimensional. The chiral condensate is, however, 
not spiral but homogeneous for small B because of contributions from all non-zero 
Landau levels. That is, the genuine chiral condensate is (y/y/) = Y.n(VnVn), among 
which only the Landau zero-mode has a special structure as in Eq. (1191 . 1 would con- 
jecture, hence, that there is no sharp phase transition from the homogeneous chiral 
condensate at B = to the chiral magnetic spiral at B ^ 0, but it may be possible 
that the inhomogeneous zero-mode contribution gradually develops, which exhibits 
a smooth crossover to the chiral spirals with increasing B. 



3 What is the chiral chemical potential? 

Equation (|8]l is very similar to the CME current (QJ, so that one might have thought 
at a first glance that Eq. (QJ emerges trivially from the insertion of y 5 in both sides of 
Eq. ®. The relation between Eqs. (QJ and dS) is not such simple, though. As a matter 
of fact, this point was a major source of confusions about the validity of Eq. ([TJ. 
One can readily extend the field-theoretical derivation of Eq. (O using Eq. (O in 
order to obtain Eq. ([TJ by introducing the axial vector fields A^, or the chiral gauge 
fields, Ar = (A^ + A 5 n)/2 and Ai = (A^ —A„)/2. Then, in the same manner as in 
the previous section, one can formulate the counterpart of Eq. (O associated with a 
vector rotation by /3 (x), that is, 



8S= / d 4 xj3(x) 



d,f (x) +N C £ ^. e W°F* v (x) Fpa (x) 



4 



(20) 



This leads to -/ - N c ^ f (q 2 f /2K 2 )e Qijk (A^ -A^)djA k = just as in Eq. (QJJ, and 
this is nothing but Eq. (QJ after the identification of Aq as /I5 (see the derivation D 
in Ref. (UJ). Although the derivation may look flawless, it triggered suspicious 
views of Eq. (QJ, which was first addressed by Toni Rebhan, Andreas Schmitt, and 
Stefan Strieker using the Sakai-Sugimoto model ll38ll . [It should be noted that the 
CME current had been exactly reproduced in the holographic models |49|.] 

Obviously, one has to deal with the chiral gauge theory with both Ar and Al to 
introduce /I5 in the above way, and it is well-known that the anomaly in the chiral 
gauge theory has a more complicated structure than that in the vector gauge theory. 
Roughly speaking, the anomaly is a consequence from the inconsistency between 
chiral invariance and gauge symmetry. In the vector gauge theory, usually, the vector 
current is strictly conserved due to adherence to gauge symmetry, and the anomaly 
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is seen in the axial vector channel only (see Eq. <j9j). In the case in the chiral gauge 
theory, however, there is no such strict demand from the theory and it should be pre- 
scription dependent how the anomaly may appear in the vector and the axial vector 
currents. Indeed we can clearly see from Eq. < [20b that the vector current is also 
anomalous. There are two representative results known as the covariant anomaly 
and the consistent anomaly, and they can coincide only when the anomaly cancel- 
lation holds, as is the case in the Standard Model. The authors of Ref. [38 1 claimed 
that the vector current should be free from the anomaly and the theory should ac- 
commodate the Bardeen counter-terms to cancel the anomalous terms in Eq. d20l ). 
Then, needless to say, the CME current is vanishing! 

This argument scared Harmen and me very much. In 2009 when Ref. [38] came 
out, Harmen was a postdoc in Frankfurt and I was also there as a visitor. Harmen's 
face is always very white, but he got even more whity, and we had a lot of discus- 
sions on Ref. [ 38 1 in Frankfurt with a fear that we might have made a big steaming 
mistake... At that time, neither Harmen nor I was 100% confident in Eq. (Q~|) (maybe 
Dima was?), and the necessity of the Bardeen counter-terms sounded plausible. This 
puzzle was one of the issues discussed in a RBRC workshop, "P- and CP-odd Ef- 
fects in Hot and Dense Matter" in May, 2010. One of the invited participants, Valery 
Rubakov, wrote a note to clarify this issue based on the discussions in the work- 
shop |39|. The essence in his argument is the following. If one introduces /I5 as the 
zeroth component of the axial gauge field, the CME current is gone indeed. How- 
ever, QCD and QED are not the chiral gauge theory. One should then introduce /X5 
in a different way as a conjugate to the Chern-Simons charge. Therefore, instead of 
adding a term flsYy^y 5 ^/ in a form of the covariant derivative in the Lagrangian, 
one should think of the Chern-Simons current which is deduced from 



2 



Ar cI I |^e Atvp %v(x)F pa (x) = ^ 



,2 



(21) 



in the QED sector. The term to be added in the Lagrangian is, 

2 

S cs = -J d A x^ 5 K°(x) = -jV c £ ^ J d A xe^ k Ai{x)djA k {x) , (22) 

from which Eq. (Q~|i immediately follows as a result of the derivative, f = 8S CS / 8Aj(x). 
One may worry about gauge invariance in the above prescription. It would be then 
more convenient to rewrite S cs in the following way after the integration by parts, 
that is manifestly gauge invariant, 

2 

S cs = j d A x9(t)N c ^ JU- E » v P° Fllv (x)F pa (x) , (23) 

where d()9(t) = /I5. In other words, we can say that /X5 is the time derivative of the 9 
angle in the QED sector, which was pointed out already in Ref. [11] and the idea of 
the charge separation driven by inhomogeneous 9 can be traced back to Ref. ||26| . 
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A subsequent question naturally arises; what happens if d(x) has not only temporal 
but also spatial dependence in general? The Chern-Simons-Maxwell theory with 
general G(x) provides us with the following modified Maxwell equations; 



v.E = P +jv c £^2(ve).B 



VxB-5 E=j+iV c ^^ [(d 0)B - (V0) x E 



(24) 



(25) 



and Faraday's law and GauB's law are not altered. We see that the CME current 
appears in the right-hand side of Eq. ( f25T > as if it is a part of the external current. 
In this manner we can conclude from Eq. (1241 that an electric-charge density is 
induced by spatially inhomogeneous (x) in the presence of B. To the best of my 
knowledge Eqs. (124-b and (f25j) are the quickest derivation of the Chiral Magnetic 
Effect, as discussed first in Ref. ||27l . 

[After I finished writing this article, I was informed by Toni, one of the authors 
of Ref. [38 1, that the confusion about the CME in the holographic context seems to 
continue. I am not able enough to make any judgment here, and those who want to 
dive into this confusion can consult the recent analysis in Ref. Q.] 



4 What really flows? 

To tell the truth, I have never gotten any satisfactory answer to the following ques- 
tion; what really flows? I have had various discussions with people who have vari- 
ous backgrounds, but those discussions ended up with more confusions than before. 
Thanks to useful conversations, nevertheless, my eyes have been open to various 
views of Eq. (fl}. People (including me) say that the CME current is an electric cur- 
rent induced by B just like Ohm's law with the electric field E. Let me begin with 
suspecting this interpretation that people just take for granted. 

In classical electrodynamics Eq. d25l ) is usually written in a slightly different way, 
i.e., 

VxB=j + a E+JV c £^ 

and d()E is called the displacement current. We see that the CME current should be 
a genuine current if <9oE can be regarded as a real electric current, for they enter 
Ampere's law on equal footing. In other words, if the displacement current is not a 
real current, the CME current is not, either. Now, we know from our experience that 
<9oE is only the time derivative of the electric field and no electric charge flows asso- 
ciated with the displacement current. The displacement current certainly plays the 
equivalent role as j as a source to create B, but it is clear that there is no movement 
of electric charge at all. It would be therefore a legitimate claim to insist that the 



(d o 0)B-(V0) xE 



(26) 
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charge separation from the CME current might be an illusion. I would emphasize 
the importance to distinguish the current and the charge in the argument here. For 
example, the most well-known example of the displacement current is the problem 
of the capacitor that is composed of two separate conductors. Let a capacitor be 
connected to the wire with finite electric current. More and more electric charge ac- 
cumulates on the conductors and produces stronger and stronger electric field inside 
as the time goes. Then, even though two conductors are physically separate and no 
electric current flows between them, the displacement current flows as if the electric 
current flowed along the wire without the capacitor. The distribution of the electric 
charge stored on the conductors is, however, totally different depending on the sit- 
uation with and without the capacitor. In this sense, thus, the charge itself may not 
flow and the charge separation may not occur with the CME current also. 

A related criticism against the CME current is that the current computed in 
Eq. (O for example is the expectation value of the current operator, y/y^xjf, and 
it is not necessarily the current. In fact, there are some studies on the Chiral Mag- 
netic Effect in the lattice gauge theory; the correlation functions of the chirality and 
the current were measured in Ref. [6], and later Eq. (HJ was checked directly on 
the lattice [48 1. It is not so straightforward, however, to interpret these lattice results 
properly. A system with a finite electric current could be steady but is out of equilib- 
rium. What one can calculate in the thermal system in equilibrium like the situation 
of the lattice simulation in Euclidean space-time is the electric-current conductivity 
according to the Kubo formula. It is a tricky question what (y/y^ y/} really represents 
in the lattice simulation. Let me take one example for concreteness. If the system 
has a condensate of the omega meson, CO* 1 , the interpolation field of co^ is ~ xjry^y/ 
and then = (yy^y) ^ 0, but this does not necessarily mean that the system has 
a persistent current. To make this point clearer, the spin operator in terms of the 
Dirac matrices is S z = ^ [y*, y y ] = jdiag[a 3 , <7 3 ], so that the spin expectation value 
is S z = (yfry) = ^(^C 3 ^) + '$r) > while the current expectation value is 

j z = (yYy) = (^<7 3 0fl) — (02 cj3 0l)' where <pL and 0# are two-component spinors 
in the left-handed and right-handed chirality, respectively. Here, the similarity be- 
tween S z and f implies that we can regard j z as a static quantity like the spin S z , 
which may well be the most appropriate interpretation of the lattice measurement. 

From the point of view of the theoretical treatment of the electric current, the 
formulation based on the linear response theory must be a good starting point. I 
believe that the work along this line in Ref. |29| should be one of the most impor- 
tant literature to think of physics of the Chiral Magnetic Effect. They computed 
the one-loop diagram on top of the ^5 background to find the chiral magnetic 
conductivity a z (co,p). The result is consistent with Eq. (T} in a particular limit; 
(J z ((0 — 0,p — ¥ 0) = liirip^oO^O,/?) — A^cL/(?/M5/2^ 2 ) which correctly repro- 
duces the CME current. In view of the result of Ref. 11291 , on the other hand, it seems 
<J X ((0 — > Q,p — 0) = 0. [This latter limit is not manifestly addressed in Ref. [29], 
but it is pointed out that the conductivity drops to one third just away from co = 0. 
It seems to be vanishing from Eq. (38) and Fig. 1 of Ref. [29].] This is a problem 
because the latter limit rather than the former one is more relevant to the real-time 
dynamics. The fact that the former limit (co — first and p = next) gives the CME 
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current ([T) suggests that the CME current should be a static quantity just as mea- 
sured in the lattice simulation and thus not a genuine electric current!? 

One may still consider that the intuitive argument leading to Eq. (O should work 
anyway. My impression is also that all above-mentioned problems are just on the 
conceptual level (though I have no idea how to reconcile them) and in practice the 
CME current flows according to Eq. (Q~|) after all. Indeed if there are almost massless 
quarks in a quark-gluon plasma and a strong B is imposed on a topological domain, 
an electric current must be induced for sure. An example of the real-time calculation 
of the CME current with not jUs but a topological domain is quite instructive in this 
sense [13|. The central innovation in Harmen's idea (as discussed in Ref. flT3l ) was 
to mimic the topological domain by putting E and B parallel to each other, with 
which £^ vpa Fj+yFpo ^ 0. Then, the particle production occurs via the Schwinger 
mechanism and the produced particles are accelerated by the fields, and the elec- 
tric current is generated. The current is time dependent and the current-generation 
rate can be analytically written down. In this setup the physical origin of the CME 
current is crystal-clear! So, if anything is fishy in physics of the Chiral Magnetic 
Effect, it should have something to do with technical defects of in equilibrium 
circumstances. 

Supposing that physics of the Chiral Magnetic Effect should be robust, let us 
admit the CME current (Q]l as it is to proceed to the next question, that is actually 
the central question in this section; what really flows? 

An intuitive explanation tells us that quarks simply flow in a quark-gluon plasma. 
It is, however, based on a classical picture, and such a picture misses quantum char- 
acter that is indispensable for phenomena related to the quantum anomaly. Look 
at the derivation of the CME current in Eq. (|7). If this derivation captures the un- 
derlying physics of the CME current, the origin of the current comes from quarks 
with infinitely large momenta. Where are such fast-moving quarks in the real quark- 
gluon plasma? They may spill out from the vacuum through quantum processes, 
but how is it possible to retrieve particles with infinite momenta? Usually the quan- 
tum anomaly involves ultraviolet regions of the momentum integration as a loop of 
virtual particles, meanwhile ultraviolet particles directly participate in the physical 
observable in the CME problem. It is very hard (at least for me) to imagine that the 
current generation in such a way really happens in a physical plasma. This deliber- 
ation brings me a further doubt about the static evaluation of the CME current. 

A natural extension of this question about the origin of the CME current is 
whether it exists in the hadronic phase and, if it does, how the current appears in 
terms of hadronic degrees of freedom. Actually this question has been something 
in mind for a long time since when we published Ref. [ 1 1 1. In the hadronic phase 
an electric current should be attributed to charged pions, but pions are insensitive to 
chirality and thus ^5 or the strong 9 angle. One possible answer would be that there 
is no CME in the hadronic phase, and if so, it would be fantastic; the CME current 
can be a signature for quark deconfinement, as implied in Ref. ATI . I had heard that 
Harmen wanted to analyze the CME using the chiral perturbation theory, though he 
never worked it out. 
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Recently I have clarified what would happen in the hadronic phase and wrote a 
paper with one of my students, Kazuya Mameda [14]. Our conclusion was a very 
natural one, and a very perplexing one at the same time. 

The CME current is unchanged even in the hadronic phase, which is very nat- 
ural since the CME current has the anomalous origin that arises from ultraviolet 
fluctuations. At low energies the anomaly should be saturated by infrared degrees 
of freedom, which is sometimes referred to as the anomaly matching. This idea is 
formulated as the Wess-Zumino-Witten action and the current should be given by 
the derivative of the total effective action with respect to the gauge field. In this way 
we found that the leading-order term in the chiral Lagrangian leads to the current, 



:e(K~id >i 7t + ~7t + id^7t-) + ... , (27) 



where E^ 1 = U^d^U, 2> = {d^U)U\ and U = e^^lf* are the standard notation 
in the chiral Lagrangian. The physical meaning of the above expression is plain 
as seen from the expansion in terms of the pion fields. It is a common form of the 
probability flow in Quantum Mechanism representing the electric current associated 
with the flow of the charged pions. 

A more non-trivial contribution comes from the Wess-Zumino-Witten part, which 
leads to the current associated with the n° domain-wall [44], i.e., 

iwzw =^L^Tff VPa (^)F pa . (28) 

This current is very similar to the CME current (Q3 and 6(x) is just replaced by 
7Z°(x)/(47t 2 f K ). Although Eq. ( 1281 is not the Chiral Magnetic Effect, it would give 
us a clue to think about the physical meaning of the CME current. Finally, the 
CME current appears from the so-called contact part of the Wess-Zumino-Witten 
action [23]; 

S P= N -L-^2 £>1Vpa J d 4 *A M (jc)(d v A p (x))d CT e(jc) , (29) 

which is just equivalent to the Chern-Simons action already discussed in Eq. d22l . 
[6 in Eq. d29l i has a different normalization by 2N{ by convention.] Naturally the 
current derived from Eq. d29l should reproduce the CME current (Q~|). This is how 
one can get the CME current in the hadronic phase and my surprise lies in the fact 
that the pion dynamics is completely decoupled from the CME current. 

Because the 7T° domain-wall looks a bit more intuitive than the mystical 9 angle, 
we shall consider a possible interpretation of the current (I28t . This is certainly a 
current, but no charged pions, n ± , are involved in the formula. Then, it is as puzzling 
in Eq. (l28l how the current can flow and what really flows. 

To answer this question, let me emphasize a very useful analogue of the Joseph- 
son current in superconductivity. The Josephson junction consists of superconduct- 
ing materials and a thin layer of insulator (S-I-S) or non-superconducting metal (S- 
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N-S) sandwiched by them. There was a big debate about whether the super-current 
can flow or not through the insulating barrier. Of course, there is no Cooper pair in- 
side of the insulator, and thereby there is nothing that takes care of the super-current. 
It should have been a natural attitude to get skeptical about such a current |4|. This 
situation, a current without current carriers, is quite reminiscent of our problem 
of the CME current or the current accompanied by the 7t° domain-wall. Everyone 
knows that the Josephson current is the experimental fact today [2 |. For the Joseph- 
son current, the coherence is the most important; in superconductor the quantum 
state is characterized by a wave-function just like a problem in Quantum Mechan- 
ics. In the QCD case, also, such a coherent state is realized by the condensation of 
fields, namely, 7t°(x) in Eq. (f28b is to be regarded as a macroscopic wave-function. 
One may then raise a question; the current of the 7t° domain-wall may be okay, 
but what about the CME current? There is no coherent field but only 9 (x) that is 
not a dynamical field but just a space-time dependent parameter! This is perfectly 
a sensible question. To answer this, I would say that 9 (x) could be promoted to the 
dynamical field without mentioning on a possibility of axion, at least in the hadronic 
phase. If the system has a pseudo-scalar (and iso-scalar) condensate such as the con- 
densate of the tj° meson (that forms rj' with a mixture with 7] 3 ), it could be mapped 
to 9(x) in the chiral Lagrangian approach. Once this mapping is noticed, there is 
no longer a big conceptual difference between the current in Eq. (f28b and the CME 
current in Eq. (Q~|i. The analogy to the Josephson current may support the reality of 
the CME current, but this argument does not tell us anything about the microscopic 
constituent of the current yet. 

Equation (|28T > means that the current can exist just with the n°(x) profile and 
the magnetic field, and then the only possible carrier of the current should be the 
quark content inside of n°. Therefore, even in the hadronic phase, I must think 
that charged quarks flow to produce the electric current. Contrary to the intuition, 
there is no inconsistency with the notion of quark confinement. Regardless of the 
presence of the flow of quarks, these flowing quarks can be still confined in a big 
wave-function of the n°(x) profile. In this way, confined quarks can flow without 
breaking confinement because of the coherent background of the meson fields. 

In reality it is next to impossible to achieve such an environment with abundant 
n° that forms a condensate to test Eq. d28l > because 7C° quickly decays into pho- 
tons via the anomalous QED process. This implies that the CME current may be 
also diminished by the photon production. Indeed, Eq. (|29T > exactly describes such 
a process of 9(x) decaying into 2y. It is interesting, besides, that one y can be pro- 
vided from B in the case with background fields. More specifically, one injected y 
and another y from B can produce a 9 (or the tj° meson), that is nothing but the 
Primakoff effect (36]. The Primakoff effect has an application as a tool to detect 
the axion [43, 37], which is understandable from the above argument once 9(x) 
is augmented as a dynamical axion field. Because physics of the Chiral Magnetic 
Effect has a connection to axion physics through 9 (x) (that was actually the very 
beginning of the path toward Ref. |28|, as I mentioned), it should be naturally moti- 
vated to think of some application of the Primakoff effect in the context of the Chi- 
ral Magnetic Effect too. Then, the reverse process of the Primakoff effect, namely, 
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7(B) + 9 — > 7, should be the most relevant to the experimental opportunity. In the 
relativistic heavy-ion collision, the profile of the magnetic field B(x) can be esti- 
mated by the simulation, and the precise measurement of 7 with subtraction of the 
background from the 71° decay is available nowadays. The unknown piece in the re- 
verse Primakoff effect is the profile of the 9 (x) distribution. Needless to say, nothing 
is more important and ambiguous than the concrete distribution of 9 (x) for any at- 
tempt to perform serious computation of the CME-related phenomena. This is why 
most of works on the Chiral Magnetic Effect address only qualitative predictions. 

I think that it must be a very interesting challenge to find a condensed-matter 
counterpart in which the Chiral Magnetic Effect may be visible and testable exper- 
imentally. This is not an unrealistic desire; axion physics can be discussed in the 
so-called topological magnetic insulator [35], and why not the Chiral Magnetic Ef- 
fect? In fact, recently, there are appearing some works one after another along this 
line. 



5 My Outlook 

Many people (including me) are still working on the theoretical aspects of the Chi- 
ral Magnetic Effect and its relatives such as the Chiral Separation Effect, the Chiral 
Magnetic Wave, etc. It is highly demanded to make some firm theoretical estimation 
about the experimental observables affected by the CME and related phenomena. To 
this end, however, one needs know the early-time dynamics even before the forma- 
tion of the quark-gluon plasma. In fact, the most interesting regime that has the 
greatest impact on the CME happens to be the most difficult regime to describe 
theoretically. 

A missing theoretical link between the coherent wave-function right after the 
collision and the thermalized plasma is the last piece of the jigsaw puzzle. Parti- 
cles should be produced from quantum fluctuations on top of coherent fields (i.e. 
the Color Glass Condensate; CGC), that translates into the entropy production. It 
is known that the coherent background fields accommodate topological flux-tubes 
that play a role of Q w in Eq. ©. Produced particles inside of those flux-tubes un- 
der a strong B should have a characteristic momentum distribution and this would 
embody the Chiral Magnetic Effect in a quantitative way. In my opinion, thus, the 
early-thermalization problem must be resolved even before talking about the obser- 
vational possibility of the Chiral Magnetic Effect, or they may have something to 
do with each other. This is because, as I stated in the previous section, the real-time 
dynamics of the Chiral Magnetic Effect should involve the Schwinger process of 
particle production, and the particle production as in the Lund string model should 
be responsible for the entropy production from fields, and thus thermalization ulti- 
mately. 

One may claim that the complete isotropization and thermal equilibrium should 
be no longer required to account for the experimental observation. This is true in- 
deed, and this is a good news for the CME, for the thermalization means that the 
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system should lose any memory and have only one information, i.e. the tempera- 
ture. If the thermalization is incomplete, it would enhance the chance that the ob- 
served distribution of particles may still remember the early-time environment like 
the presence of the strong B and/or the topological flux-tubes. For the purpose of 
testing the idea as compared to the experimental data, it is indispensable to perform 
some serious simulation of the early-time dynamics of the heavy-ion collisions. 

Unfortunately, there is no successful simulation starting from the CGC initial 
condition to achieve some reasonable input for the hydrodynamics within a reason- 
able time scale (See Ref. [17| for a latest attempt). There are so many theoretical 
efforts in this direction including mine [ 1 1 and it should be definitely worth dis- 
cussing them, but not here and on another occasion maybe. In this article I have 
discussed physics of the Chiral Magnetic Effect and presented my views on the 
physical interpretation. Now my story has become a bit too diverging, and I should 
stop here, with a hope that some readers may find my views useful for future inves- 
tigations. 
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